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The adiabatic projection method is a general framework for studying scattering 
and reactions on the lattice. It provides a low-energy effective theory for clusters 
which becomes exact in the limit of large Euclidean projection time. Previous stud¬ 
ies have used the adiabatic projection method to extract scattering phase shifts from 
finite periodic-box energy levels using Liischer’s method. In this paper we demon¬ 
strate that scattering observables can be computed directly from asymptotic cluster 
wave functions. For a variety of examples in one and three spatial dimensions, we 
extract elastic phase shifts from asymptotic cluster standing waves corresponding to 
spherical wall boundary conditions. We find that this approach of extracting scat¬ 
tering wave functions from the adiabatic Hamiltonian to be less sensitive to small 
stochastic and systematic errors as compared with using periodic-box energy levels. 

PACS numbers: 21.60.De,04.60.Nc,25.55.-e 


I. INTRODUCTION 


Ab initio description of scattering and reactions involving nuclei is one of the major chal¬ 
lenges in computational nuclear physics. Recent progress along this line has been achieved 
using resonating group methods hhbi, fermionic molecular dynamics m is], the coupled- 
cluster expansion, see [6] for a review article, and variational and Green’s function Monte 
Carlo methods [7, 8J. For calculations involving lattice methods, there has been progress in 
using finite periodic volumes to analyze coupled-channel scattering [SHE] and three-body 
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systems [IMS]. In this paper we concentrate on the adiabatic projection method for cal¬ 
culating nuclear reactions from lattice simulations in the framework of chiral effective field 
theory (EFT). See [HH22] for some recent results using lattice EFT. The general strategy 
in the adiabatic projection formalism as formulated in the pioneering work [23H25j involves 
two steps. First, one uses the Euclidean time projection method to determine an adiabatic 
Hamiltonian for the participating nuclei starting from the microscopic Hamiltonian derived 
in chiral EFT. In the second step, one uses a technique such as Luscher’s finite-volume 
method [261127] to extract the corresponding scattering phase shifts [2.4] . 

The adiabatic projection formalism has been successfully benchmarked against continuum 
calculations for fermion-dimer scattering. However for heavier systems, Luscher’s finite- 
volume energy approach for extracting scattering phase shifts is expected to suffer from 
potentially large errors due to stochastic and systematic uncertainties in the lattice Monte 
Carlo energies. In this paper, we explore various techniques to access scattering on the 
lattice that do not require a high-accuracy determination of the energy spectrum. For 
simple three-body systems in one and three dimensions, we demonstrate that scattering 
phase shifts can be reliably extracted from the asymptotic cluster wave functions within the 
adiabatic projection method. 


Our paper is organized as follows. We begin with introducing the adiabatic projection 
method in some detail in section [H} Section 0 describes Luscher’s finite volume method, 
while section [IV] describes the extraction of the asymptotic cluster wave functions on the 
lattice. The microscopic Hamiltonian used in our work is specified in section [Vj which also 
provides details on the computation of the adiabatic Hamiltonian. Various approaches for 
extracting the two-cluster elastic scattering phase shifts in one and three spatial dimensions 
on the lattice are introduced and applied in sections [Vi] and |VII[ respectively. Finally, the 


main results of our study are summarized in section VIII 


II. THE ADIABATIC PROJECTION METHOD 


The adiabatic projection method treats the cluster-cluster scattering problem on the lattice 
by using Euclidean time projection to determine an adiabatic Hamiltonian for the partici¬ 
pating clusters. When the temporal lattice spacing is nonzero, an adiabatic transfer matrix 
rather than the Hamiltonian is constructed, but the method is essentially the same. We 
start with an L 3 periodic lattice and set of two-cluster states | R) labeled by their separation 
vector R , as illustrated in Fig.[l In general, there are spin and flavor indices for these states, 
but we suppress writing the indices for notational simplicity. The exact form of these two- 
cluster states is not important except that they are localized so that for large separations 
they factorize as a tensor product of two individual clusters, 

ifl> = X>+(i) 

r 

These states are propagated in Euclidean time to form dressed cluster states, 

\R) T = exp(-Hr)\R). 


( 2 ) 
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Figure 1: (Color online) Two-body cluster initial state | R) separated by the displacement vector 
R. 

By evolving in Euclidean time with the microscopic Hamiltonian, deformations and polariza¬ 
tions of the interacting clusters are incorporated automatically, and we are projecting onto 
the space of low-energy scattering states in our finite volume. In the limit of large Euclidean 
time, these dressed cluster states span the low-energy subspace of two-cluster continuum 
states. 

We evaluate matrix elements of the microscopic Hamiltonian with respect to the dressed 
cluster states, 

[ffr]** = r(R\H\S) r . (3) 

Since the dressed cluster states \R) T are, in general, not orthogonal, we construct a norm 
matrix N T given by the inner product 

[Wr]** = r(R\R')r- (4) 

From this a Hermitian adiabatic Hamiltonian matrix can be defined using the inverse square 
root of the norm matrix, 

[»;]«,«"= E [ JV r- 1/2 ] 4S »[«r]*,,*dJV- ,/2 ] S „, a ,. (5) 

R",R!" 

In the limit of large r, the spectrum of exactly reproduces the low-energy finite volume 
spectrum of the microscropic Hamiltonian H. So for the elastic phase shifts, one can take the 
spectrum of and apply the finite-volume scaling analysis developed by Liischer [261 EH- 


III. LUSCHER’S FINITE-VOLUME METHOD 


The Liischer method [261 EH relates the two-body scattering states in periodic finite volume 
to the scattering parameters in the infinite volume and continuum limits. The main idea 
behind the derivation of this method is to use the known asymptotic form of the wave 
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function for a short range interaction and the periodicity of the system. The one-dimensional 
result of this ansatz is [ 2 f>] 


e 2i<5(p) _ e ~ipL 


( 6 ) 


where p is the relative momenta between two clusters 1 , 5(p) is the scattering phase shift, and 
we assume the total momentum of the two-cluster system to be zero. In three dimensions, 
the situation is more complicated due to breaking of the rotational invariance by the cubic 
symmetry of finite periodic box. The scattering phase shifts are directly related to the 
momentum via the formula [2EI, 27j 


p cot 6e(p) = — yS(rj) for £ = 0 , 1 , 

7T L 


where rj = ) 2 and S{ij) is three-dimensional zeta function, 


S(rt) = lim 

A->oo 


^ 0(V-^) _ 4irA 


n — tj 


(7) 


( 8 ) 


or, in the exponentially accelerated form [I26H28] . 


S(p) = 2 tt 3 / V(2?7 - 1) \ c" 


3 -|n | 2 


-7T 3 / 2 


\n | z — r / 


r 1 e Xr i 

i dX W 2 


4A 2 ?7 2 — e _7r2 l"l 2 / A 


(9) 


The relation between the relative momentum appearing in Eq. 0 -@ and the finite-volume 
energies for £ = 0 is given by [35j 36j 


2 

E(p, L) = ^--B 1 -B 2 + rxO)A E^L) + f 2 ( V )AE 2 (L) , (10) 

where p is the reduced mass of the system, B, is the binding energy of the cluster i = {1,2} 
in the infinite volume limit, A Ei(L) = E t (L) + Bi is the finite volume energy shifts of the 
clusters in the rest frame, and fi(rj) is the topological correction factor to the energy of the 
cluster i, 


r(p) 


Xfc {k 2 - v) 


-2 


E 


XL cos qQ 

3 (k 2 — r/) 2 


( 11 ) 


In Ref. 


it was found that for £ > 0 , the topological corrections are suppressed by the 


size of the finite volume, r(? 7 ) = 1 + 0(1/L), so that Eq. (10) becomes 


E(p,L) 


2/i 


+ Et(L) + E 2 (L). 


( 12 ) 


1 Clusters refer to either a point-like particle or a composite particle as a bound state of several particles. 
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Figure 2: (Color online) The colored data points show the lowest energy levels of the adiabatic 
Hamiltonian versus projection time r for the dimer-fermion system on various L 3 periodic lattices. 
The horizontal black-dashed lines are the lowest energy levels of the microscopic Hamiltonian for 
the same L 3 lattices. The dimer energy is —2.2246 MeV and the lattice spacing is 1.97 fm. 


To minimize potential errors in the calculations using the Liischer’s method, we also take 
into account the effective mass of the clusters on the lattice. Luscher’s method is a useful and 
commonly used tool for calculating scattering parameters. For examples of recent extensions 
and generalizations of Liischer’s original work see Ref. HEBEI EEH2]. 

The simple elegance of Liischer’s method is that all of the information regarding scattering 
phase shifts is encoded into finite-volume energy values. This simplicity can, however, be a 
weakness when applied to scattering processes relevant to low-energy nuclear physics. The 
problem is that the binding energies of the scattering nuclei can be anywhere from a few 
MeV to tens or hundreds of MeV, while the finite-volume scattering energy shifts can be as 
small as a few keV. The problem is even more difficult in lattice QCD calculations where 
the rest energy of the nucleons is also part of the calculations. 

Thus, while the scattering data is encoded in the finite-volume energy and waiting to be 
extracted, the finite-volume energy value is prone to several sources of potentially large 
errors. In Fig. [2| we show the lowest-lying energy state for a dimer-fermion system versus 
projection time r for various L 3 periodic lattices. The dimer energy is set to be —2.2246 MeV 
and the lattice spacing is 1.97 fm. We consider this example in detail later in our discussion. 
The point we emphasize here is that in order to measure the s-wave scattering phase shift 
to an error of a few degrees, we need to measure the finite-volume energy to an accuracy of 
about 10 keV. In this simple three-particle calculation we are using exact matrix methods, 
and there are no stochastic errors. However, in a typical large-scale Monte Carlo calculation, 
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there are stochastic errors that grow exponentially with projection time r. The stochastic 
errors are much reduced in a constrained Monte Carlo calculation using, for example, fixed 
fermionic nodal constraints. However here we have a different problem that the scattering 
energies may be artificially shifted by the constraints. In addition to these issues, there 
are also corrections to the binding energies of the scattering nuclei due to the finite volume 

mm- 

In view of the problems with finite-volume energy calculations for low-energy nuclear scatter¬ 
ing, we introduce in this paper another approach for extracting scattering phase shifts which 
directly analyzes cluster wave functions generated by the adiabatic Hamiltonian. This ap¬ 
proach has the advantage of being far less sensitive to small errors in reproducing the binding 
energy and detailed structure of the participating nuclei. 


IV. ASYMPTOTIC CLUSTER WAVE FUNCTIONS 


For each particle of mass m, the microscopic Hamiltonian contains a kinetic energy term of 
the form —V 2 /(2m) or, more precisely, the lattice approximation to this operator. Therefore 
the Euclidean-time evolution operator exp(— Hr) acts as a diffusion operator with diffusion 
constant inversely proportional to the particle mass m. Since the particles are diffusing in 
space as a function of Euclidean time, it is necessary to be precise about what we mean by 
widely-separated clusters at asymptotically large distances. 

In order to explain the various time and length scales of our asymptotic wave function 
analysis, it is useful to first specify a relative error tolerance, e, for all steps in our cluster- 
cluster scattering calculation. As we project to large Euclidean time, any isolated single¬ 
cluster initial state will simply relax into the ground state of that cluster system. We define 
r e as the time at which the relative contamination due to excited cluster states is less than 
e. 

We now consider applying Euclidean time projection for the time duration r e in order to 
remove excited cluster states. During the time interval r € , each cluster undergoes spatial 
diffusion by an average distance proportional to y/r e /M, where M is the mass of the cluster. 
We call this distance the diffusion length. Let d £: i be the diffusion length for the first cluster, 
and d f 9 be the diffusion length for the second cluster. In order to have some widely-separated 
clusters, we take our periodic box length L to be much larger than d cA and d Cj 2 . We recall 
that R is the initial separation vector between the two clusters. When \R\ d cA , d C 2 
the dressed cluster state \R) Te consists of non-overlapping clusters. We thus can define an 
asymptotic radius R e as the radius such that for \R\ > R, the amount of overlap between 
the cluster wave packets is less than e. 

In the asymptotic region \R\ > R f , our dressed clusters are widely separated and interact 
only through long range forces such as the Coulomb interaction. Therefore, we can describe 
our system in terms of an effective cluster Hamiltonian H eS that is simply a free lattice 
Hamiltonian for two point particles accompanied by the long-range interactions inherited 
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from the microscopic Hamiltonian. So in the asymptotic region we have 



(13) 


(14) 


where the coefficient c is given by the overlap of the initial single-cluster states in Eq. ([Tj) 
with the exact single-cluster energy eigenstates. Since 



we conclude that the adiabatic Hamiltonian coincides with the effective cluster Hamiltonian 
in the asymptotic region, 



This last result is quite significant. We are using Euclidean time projection to calculate 
the matrices N T and H r . This Euclidean time projection entails a considerable amount of 
diffusion of the clusters. However, in the asymptotic region, we are in essence inverting the 
diffusion process when computing the adiabatic Hamiitonian and are left with an effective 
cluster Hamiltonian in a position space basis. 

For cases where there are no long range interactions, the scattering states of the adiabatic 
Hamiltonian are given by a superposition of Bessel functions in the asymptotic region. For 
the case with Coulomb interactions, the scattering states of the adiabatic Hamiltonian in 
the asymptotic region correspond to a superposition of Coulomb wave functions. In the rest 
of this paper we extract cluster-cluster phase shifts for a variety of different examples in 
one and three spatial dimensions. We use spherical hard wall boundaries placed at some 
large wall radius R wa ,n in the asymptotic region and determine phase shifts by analyzing 
the asymptotic scattering wave functions. As we will show in the following discussion, this 
method is robust and accurate even in cases where the corresponding finite-volume energies 
produce large errors using Liischer’s method. 


V. LATTICE FORMALISM 


A. Microscopic Hamiltonian 


Throughout this work, we consider several systems, all of which are comprised of three point 
particles. For some cases the particles are distinguishable and in some cases not. We specify 
the quantum statistics explicitly at a later stage. We consider the limit of where the range 
of the interaction is negligible compared to the scattering length. Then, in the low-energy 
limit, the interaction between particles can be represented by a delta-function. 

Let b s and b\ be the annihilation and creation operators for each particle species s, and let 
p s (r) be corresponding density operator, 


Pa(f) = b\{f)b s (r ). 


( 17 ) 








(18) 


The Hamiltonian in D spatial dimensions has the form H = Ho + V where 

tfo = -£ f d D rbl(r)^-b s (r ), 

S 

V = ^2 C s,s’ f d D r : p s (f) p s '{r) :, (19) 

s<s' 


and the :: symbols indicate normal ordering. On the lattice, we can write the Hamiltonian 
as H = H 0 + V where 


= 2 2hs ^> ~ bs + 0 ~ b s{n~ 1) 


l n 


V = J2J2 Cs ’ s ' : d. s '(«) : • 


s<s' ft 

Here l denotes lattice unit vectors along all possible spatial axes. 


( 20 ) 

( 21 ) 


For D > 1, a regularization of ultraviolet divergences due to the zero-range interaction is 
needed and provided by the nonzero lattice spacing. We denote the spatial lattice spacing as 
a. Throughout this paper, we will write all quantities in lattice units, which are physical units 
multiplied by the corresponding power of a in order to render the combination dimensionless. 


In this paper we consider two examples, the first in one dimension and the second in three 
dimensions. We have organized our discussion to discuss both examples together, illustrat¬ 
ing the methods and results for both examples, one after the other. For both cases the 
lattice spacing is taken to be 1.97 fm or 0.0100 MeV -1 . For the one-dimensional example, 
we consider three distinguishable particles with equal masses and calculate the scattering 
between a particle of species 3 with a dimer composed of particle types 1 and 2. For this 
case we tune the coupling C \ t 2 to produce a bound state with energy —2.0000 MeV. The 
couplings Ci t 3 and G' 2.3 are both set equal to exactly one-tenth of the coupling CV- 2 - 

For the three-dimensional example, we consider two species of fermions which we label 
with spins, f an d V In this case we calculate scattering between a j" particle and dimer 
composed of particles j" and 4" The fermion-dimer system corresponds to the neutron- 
deuteron scattering in the spin-quartet channel at leading order of pionless effective field 
theory. Therefore, we set the coupling constant Cyj, to a value for which the dimer energy 
has the value of the physical deuteron energy of —2.2246 MeV. 


B. Adiabatic Hamiltonian 


We now discuss the implementation of the adiabatic projection method using a set of initial 
cluster states as introduced in Eq. ([!]). The two examples we consider are a particle-dimer 
system in one dimension and a fermion-dimer system in three dimensions. For the one¬ 
dimensional particle-dimer system, the particle-dimer cluster states are defined as 

\4>g) = V b l l . i! -) b Un) l >3(n + R) 10} ■ 


( 22 ) 
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Figure 3: (Color online) Left panel: Finite-volume energies extracted from the microscopic Hamil¬ 
tonian and the two-cluster adiabatic Hamiltonian for the particle-dimer system in one dimension. 
Right panel: The particle-dimer scattering phase shifts calculated from the data in the left panel 
using Liischer’s method. Color online: Red (blue) symbols show the results corresponding to the 
original (adiabatic) Hamiltonian. 

We use vector notation for notational consistency even though there is only one spatial 
direction. For the three-dimensional fermion-dimer system, we write the fermion-dimer 
initial cluster states as 


I (Pr) = b\(n)b\(n)b\(n + R) |0), (23) 

n 

where R ^ 0 because of Fermi statistics. 

We now project the initial cluster states \(f>g.) in Euclidean time r using the microscopic 
Hamiltonian, H , for the respective systems. We use the Trotter approximation, 

exp(—Hr) « (1 - a t H ) Lt , (24) 

where a t is a time step parameter and r = a t L t . For both the one-dimensional particle-dimer 
and three-dimensional fermion-dimer systems we have 

10g> T = (1 - a,H ) L * |^>. (25) 

For large r we obtain an accurate representation of the low-energy spectrum of H using the 
adiabatic Hamiltonian defined in Eq. (|5j). 


C. Scattering phase shifts from periodic-volume energy levels 


In Fig. [3| we compare the energy spectrum of the microscopic Hamiltonian and the two- 
cluster adiabatic Hamiltonian for the particle-dimer system in one dimension for r = 
0.30 MeV -1 . We use the hnite-volume energies to calculate the particle-dimer scattering 
phase shifts employing Luscher’s method. Comparative results for the scattering phase 
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Table I: The lowest-lying energy and the corresponding scattering phase shifts for the dimer- 
fermion system computed using the Luscher’s method for various periodic lattices. The energies 
E correspond to the microscopic Hamiltonian while the energies E(r) correspond to the adiabatic 
Hamiltonian with projection time r. 


L 

H 




E (MeV) 

5 (degrees) 

r (MeV 1 ) 

E(t) (MeV) 

5 (degrees) 

8 

-1.4423319 

-42.6 

0.37 

-1.4060289 

-44.1 

9 

-1.6670941 

-37.9 

0.37 

-1.6121233 

-40.9 

10 

-1.8171997 

-33.6 

0.34 

-1.7214154 

-40.8 

11 

-1.9203247 

-29.8 

0.34 

-1.8054714 

-39.6 

12 

-1.9929256 

-26.4 

0.34 

-1.8617182 

-40.0 


shifts are shown in the right panel of Fig. [3j Although the energy spectra of the adiabatic 
Hamiltonian and microscopic Hamiltonian are very similar, the resulting phase shifts have 
large differences at low energies. We also see a disagreement at low energies among phase 
shifts determined using different adiabatic Hamiltonian energy levels. This is because low- 
energy phase shifts are computed using very large box sizes L where the level spacing is 
small, and this magnifies any small discrepancies in the energy values. 

We have performed a similar analysis for the fermion-dimer system in three dimensions. In 
Table [TJ we compare the lowest-lying energy states and resulting phase shifts computed using 
the microscopic Hamiltonian and adiabatic Hamiltonian for the dimer-fermion system. The 
energies E correspond to the microscopic Hamiltonian, while the energies E{t ) correspond 
to the adiabatic Hamiltonian. The scattering phase shifts are calculated using Luscher’s 
method. As can be seen from Table [IJ the accuracy of the finite-volume energies is the range 
of 40-130 keV. However the relative error in the resulting phase shifts becomes as large as 
50%. 


Fortunately the adiabatic Hamiltonian contains more usable information than just periodic- 
lattice energy levels. As we show in Section [VT| the scattering phase shifts can be determined 
with far better accuracy using the properties of the asymptotic scattering wave function. 


VI. SCATTERING CLUSTER WAVE FUNCTION: METHODS 


Borasoy et al. introduced a method to compute phase shifts for point-like particles on a 
lattice using a spherical wall boundary I32J- As done in the continuum |3_8j , a spherical 
hard wall of radius i? wa ii is imposed on the relative separation of the two particles. In this 
study, we consider two-cluster systems, and the spherical hard wall boundary is imposed on 
the relative separation of the two clusters. For two clusters interacting via a potential of a 
finite-range R, the wave function at distances r > R is given by 

= An cos (pr + 5g — Ur/2) for one dimension, 

= R^\r) Yn :Tni ,(6, 0) for three dimensions, 


(26) 

(27) 
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Table II: Irreducible representations of the cubic rotation group SO(3, Z) and relation to spherical 
harmonics for £ < 2. 


SO(3, Z) 

SO(3) 


At 

e = o 

{%>} 

T\ 

£ = \ 

We y M } 

Et 

1 = 2 

{Y 2 ,0 , (W,-2 + Y 2 , 2 )/V 2} 

t 2 

£ = 2 

{Y 2 , i, (W,-2 - Y 2 , 2 )/V2, Y 2 -i} 


where p is the relative momentum of the clusters. For the one-dimensional case, there is no 
angular momentum, but we nevertheless use the notation £ = 0 for even parity and £ = 1 
for odd parity. In three dimensions, the total wave function is decomposed into the radial 
part R^\r) and spherical harmonics Y^ me {9 ,0). The radial wave function R^\r) has the 
asymptotic form 


R { [\r) = A t [cos Se(p) je(pr ) - sin 8tip) n e (pr)] , (28) 

where At is a normalization coefficient, and jt and rq denote sphericai Bessel functions of 
the first and second kinds. Therefore, the three dimensional wave function in Eq. (27) can 
be rewritten as 


(r) = AtY^ mi (0, 0) [cos 5tip) jt(pr) - sin 5tip) n e (pr )] . (29) 


In the asymptotic region, we fit Eq. (29) to the lattice wave functions emerging from imposing 


the spherical hard wall at radius -R wa ii. We evaluate the spherical harmonics Y( jnt {9, 0) on 
the lattice points, noting that there is no exact separation of radial and angular variables on 
the lattice. One must take into account the break up of the 2£ + 1 spin multiplets according 
to irreducible representations of the cubic rotation group [HH1 SDj • Irreducible representations 
of the SO(3, Z) cubic rotation group are given in Table [TT] for £ < 2 [JlJ. 

For £ = 0, the spherical harmonic To i0 (d, 0) is angle-independent, and we can directly match 
Eq. (29) to the lattice wave functions. However, for £ > 0, the angular dependence makes 


the fitting more difficult. To resolve this issue, we use the identity 


E IW = 

rri£=—£ 


2 2 £ +1 


47r 


(30) 


Since the angular dependence drops out of this expression, it is convenient to work with the 
wave function probability distribution summed over rrq, 



2 £+ 1 

47r 


\A e \ 2 [cos S e ip)jtipr) 


sin 5 e (p) n e (pr)] 2 . 


(31) 


To obtain accurate results for the phase shifts within this approach, it is helpful to address 
some uncertainties in the precise location of the wall radius R wa \\- While we impose a very 
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r (lattice units) 


Figure 4: (Color online) First three lowest-lying s-state wave functions (left panel) and p -wave 
probability distributions (right panel) for the fermion-dimer system in three dimensions. A spher¬ 
ical hard wall is imposed at 2? wa ii = 13.0 lattice units. 


large but finite repulsive potential at distances r > -R wa ii, a close examination shows that 
wavefunction vanishes at some slightly larger radius /?' vall = i? wa n + e, where e is some 
fraction of a lattice spacing. This is illustrated in Fig. [4] for the lowest three s- and /estates 
of the fermion-dimer system in three dimensions. The hard wall potential is imposed for 
r > R wa ,\\ = 13.0 lattice units. We find that the first zero of the Is wave function is at 13.32, 
the second zero of the 2s wave function is at 13.32, and the third zero of the 3s wave function 
is at 13.31 lattice units. For the first three lowest /estates, we find that corresponding zeros 
are all at 13.32 lattice units. 

One can extract the scattering phase shifts by doing a three parameter fit of the overall 
normalization, momentum and phase shift to the interacting wave functions. We present 
the results of this fitting procedure later in our discussion. However, we have found more 
accurate results by making use of the empirical observation that R' wai] changes very little 
when going from the non-interacting system to the interacting system at approximately the 
same scattering energy. We first determine /?( vall from the lattice wave functions of the non¬ 
interacting cluster-cluster system. Then, using the same value of /?( vall , we fit the interacting 
wave functions using a two parameter fit to determine the phase shift of the interacting 
system using the relations 


St(p) 


-PK 


'wall 


tt(1+1) 


mod 7T 


tan 


-l 


MpR wall/a) 


for one dimension 
for three dimensions. 


(32) 


In Fig.[5j we show the Is and 2s non-interacting particle-dimer wave functions used to calcu¬ 
late R' wai] and the corresponding interacting fermion-dimer wave functions used to determine 
the s-wave scattering phase shift 5 0 (p). 
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r (lattice units) r (lattice units) 


Figure 5: (Color online) Matching of the wave functions of the Is- (left panel) and 2s-state (right 
panel). Free and interacting wave functions are denoted by (color online: blue) saltires and (color 
online: red) crosses, respectively. 


VII. SCATTERING CLUSTER WAVE FUNCTION: PHASE SHIFT RESULTS 


We now compute the scattering phase shifts using the adiabatic projection method and 
scattering cluster wave functions in our one-dimensional particle-dimer system and three- 
dimensional fermion-dimer system. The results are benchmarked against phase shifts ex¬ 
tracted from the exact three-body energy spectrum obtained using Liischer’s method. For 
the three-dimensional fermion-dimer system, the three-body energies are computed using 
the Lanczos iterative eigenvector method with a space of L 6 basis states. These can be 
viewed as exact lattice phase shifts. We note that while the adiabatic projection method 
calculations can be applied to much larger systems using lattice Monte Carlo, these exact 
Lanczos calculations are limited to small systems. 


A. Particle-dimer scattering in one dimension 


We make use of the simplicity of the one dimensional system to benchmark several different 
techniques for obtaining phase shifts. In each case we construct the adiabatic Hamiltonian 
using a projection time of r = 0.30 MeV. Luscher’s energy spectrum method makes use of 
the energy of a scattering state and not its wave function. However we can also calculate 
phase shifts by fitting the wave function in the periodic box to its asymptotic form Eq. (26). 
Here and in what follows, we refer to this method as the Luscher wave function method. 
In Fig. [6|a) we show an example of such a fit. The resulting value of the phase shift 
in this example is S 0 (p) = —160.5 ± 0.3° with momentum p = 31.28 ± 0.03 MeV. We can 
compare the phase shifts obtained with this approach to Luscher’s method using the exact 
energy spectrum. Fig. |6](b) shows Luscher’s energy spectrum method phase shifts obtained 
for a number of lattices with L = 6... 100 in lattice units. With the Luscher periodic- 
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Figure 6: (Color online) Particle-dimer phase shifts in one dimension calculated using the Liischer 
wave function method. Left panel: An example of the wave function matching. Right panel: 
Comparison of the phase shifts calculated using the Liischer periodic-box wave function method 
and Luscher’s finite-volume method with the exact energy spectrum. 


box wave function method, phase shifts are calculated for L = 10, 20,... 100 in lattice 
units. As expected, one observes very good agreement between the Liischer energy spectrum 
method and the Liischer wave function fit. While this method clearly works very well for 
the one-dimensional system, we find that fitting periodic-box wave functions is much more 
problematic and less accurate in three dimensions, especially for £ > 1. This is most likely 
due to systematic errors arising from lattice spacing artifacts. 

The second method we consider is the spherical wall approach described in the previous 
section. We impose a hard boundary on the relative separation and calculate the phase shifts 
from the properties of the standing wave functions. This method allows for the calculation of 
several data points per chosen lattice volume, since we can vary the value of the wall radius, 
_R wa ii. This represents an important computational advantage as compared to the Liischer 
wave function method, especially for calculations in three dimensions. We will consider two 
versions of the spherical wall method. In the first version we do a three parameter fit of the 
overall normalization, momentum and phase shift of the interacting wave functions. Fig. [7] 
shows the application of this approach for the case L = 50 and i? wa ii = 20 in lattice units. The 
resulting value of the phase shift in this example is So(p) = —163.0 ± 0.4 for the momentum 
p = 36.17 ± 0.07 MeV. This method also shows very good agreement with Luscher’s energy 
spectrum method. The phase shifts are calculated for L = 50 and -R wa ii = 13 ... 23 in lattice 
units. 

Our next approach is a second version of the spherical wall wave function method. In this 
case we determine R( vall = R w a ii + e from the non-interacting particle-dimer wave function. 
In the example shown in Fig. [8j the boundary is set at R Wil \\ = 17, and we find the wave 
function vanishes at i2^ all = 17.901, and the momentum of the free wave function is po = 
43.874 ± 0.02 MeV. We then do a two-parameter fit to the interacting wave function and 
find So(p) = —165.8 ± 0.5 for the momentum p = 42.5 ±0.1 MeV. The phase shifts shown in 
the right panel of Fig. [8] are calculated for L = 50 and R wa n = 13 ... 23 in lattice units. The 
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Figure 7: (Color online) Particle-dimer phase shifts in one dimension calculated using the spherical- 
wall method in one dimension with a three-parameter fit. Left panel: An example of the wave 
function fits. Right panel: Comparison of the phase shifts calculated using the spherical wall 
method and Liischer’s finite-volume method with the exact energy spectrum. 
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Figure 8: (Color online) Particle-dimer phase shifts in one dimension using the second spherical wall 
approach with R' wdb n determined from the non-interacting wave function. Left panel: An example 
of the wave function fits. Right panel: Comparison of the phase shifts calculated using the second 
spherical wall approach and Liischer’s finite-volume method with the exact energy spectrum. 


results are in agreement with the results of the first spherical wall approach, but have smaller 
error bars, especially for the odd parity phase shifts. In three dimensions this improvement 
becomes more significant. 

One of the disadvantages of the spherical wall method is that one needs to go to rather 
large values of R wa \\ and L in order to probe very low energies. The last method we consider 
overcomes this issue. In order to compute phase shifts at low momenta using small lattices, 
we impose a spherical hard wall and add also an attractive well potential in front of the 
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Figure 9: (Color online) Particle-dimer phase shifts in one dimension with the spherical wall and 
an attractive well potential. Left panel: An example of the wave function matching. Right panel: 
Comparison of the phase shifts calculated using the approach based on the combination of the 
spherical wall and attractive well with Liischer’s finite-volume method using the exact energy 
spectrum. 


wall boundary. We treat the depth of the well as an adjustable continuous parameter. The 
example shown in Fig. [9] corresponds to the case of L = 30, i? wa ii = 13 and R we \\ = 12 in 
lattice units. The resulting value of the phase shift in this example is S 0 (p) = —155.0 ± 1.6° 
at momentum p = 24.9 ± 0.6 MeV. The phase shifts shown in the right panel of Fig. [9] 
are calculated for L — 50, i? wa ii = 23 and different well depths. The agreement with 
Liischer’s energy spectrum method is very good, and we also obtain phase shifts for smaller 
momenta. However, the additional attractive potential distorts the asymptotic form of the 
wave function near the potential well. The distortion of the wave function grows with the 
depth of the attractive potential, and the calculation of the phase shifts for smaller momenta 
is achieved at the expense of a larger relative error. This might complicate the application 
of this method for calculations in three dimensions where the values of R W ai\ are typically 
smaller. 


B. Fermion-dimer scattering in three dimensions 


Before presenting lattice calculations for fermion-dimer scattering in three dimensions, we 
first review continuum calculations of the same system in the limit of zero range interactions. 
This corresponds to neutron-deuteron scattering at leading order in pionless effective field 
theory [42VI44] . The T-matrix is obtained by solving the Skorniakov-Ter-Martirosian (STM) 
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Figure 10: (Color online) The s-wave (left panel) and p-wave (right panel) scattering phase shifts 
for fermion-dimer scattering in three dimensions. We compare phase shifts calculated using the 
the spherical wall approach with R' wdb n determined from the non-interacting wave function with 
Luscher’s finite-volume method using the exact energy spectrum. 


integral equation, 

T e (k,p) = -^-Qe( 
mpk \ 

2 r°° q TAk, q ) _ f p 2 + q 2 — mE — «0 + 

-/ dq - - Qe - 

7T Jo P \/3q 2 /4: - mE - *0+ - 7 V PQ 


p 2 + k 2 — mE — z0 + 
pk 


(33) 


where 7 is the dimer binding momentum, E = 3 p 2 / (4m) — 7 2 /m is the total energy, and Qe 
is the Legendre function of the second kind, 


Qe{a) 


1 

2 


r dx Ep. 

J—l X + CL 


(34) 


The scattering phase shift can then be extracted by matching the solution of Eq. (33) with 
the on-shell T-matrix 


Tt(p,p) = 


37T 


P‘ 


2 1 


m p 2£+1 cot §£ — ip 2£+1 


(35) 


For the three-dimensional fermion-dimer calculation, we use only the most promising of the 
four approaches explored previously for the one dimensional system. This is the spherical 
wall method with R ' wau determined from the non-interacting wave function. In Fig. 10, the 


squares show lattice results for the s-wave and p-wave scattering phase shifts using adiabatic 
projection method with r = 0.37 and spherical wall method. The circles are the exact 
lattice results obtained using Luscher’s method applied to the energies of the microscopic 
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Hamiltonian H evaluated using Lanczos eigenvector iteration. The dashed lines correspond 
to leading order pionless EFT continuum results obtained using the STM equation. The 
solid lines are fits of the lattice data using an effective range expansion, 

p 2e+1 cot 8(.{p) = + \r t p 2 + 0{p 4 ). (36) 

ae 2, 

We note again that these exact Lanczos benchmark calculations using the energies of H are 
only possible in small systems. The adiabatic projection method is needed to probe much 
larger systems via lattice Monte Carlo. We have seen in Table [I] that when we use Luscher’s 
method to extract phase shifts from the adiabatic Hamiltonian energies, the errors are as 
large as 50% at low energies. In comparison with this, we observe much smaller error bars 
and excellent agreement between the adiabatic Hamiltonian results and the exact lattice 
phase shifts in Fig. [lOj We should mention that the discrepancies between the lattice and 
continuum results at large momenta are nothing more than lattice spacing artifacts and 
would go away in the limit of vanishing lattice spacing. 


VIII. SUMMARY AND CONCLUSIONS 


In this paper we have presented a new method for computing scattering parameters directly 
from cluster wave functions. It allows us to bypass computation of the energy spectrum 
and thus to avoid potentially large errors in calculating low-energy nuclear scattering and 
reactions using the adiabatic projection method. We showed that the adiabatic Hamiltonian 
in the asymptotic region reduces to a simple cluster Hamiltonian in a position space basis. 
In this way we can extract scattering phase shifts directly from the scattering cluster wave 
functions. 


We considered particle-dimer scattering in one dimension and fermion-dimer scattering in 
three dimensions. In the one dimensional particle-dimer example, we explored various ver¬ 
sions of the adiabatic projection method using the cluster wave functions rather than the 
finite volume energies to extract phase shifts. First, we presented a simple matching of 
the adiabatic lattice wave function in periodic box to the asymptotic form in Eq. (26). 
We were able to accurately fit the wave function in the asymptotic region, and the calcu¬ 
lated phase shifts were in good agreement with exact lattice phase shifts computed using 
Luscher’s energy method and the spectrum of the microscopic Hamiltonian. Next, we im¬ 
posed a spherical hard wall boundary on the relative separation and matched the resulting 
wave function to Eq. ( p6| ) . This approach was also found to agree accurately with the exact 
lattice phase shifts. We considered two variants of this spherical wall approach and found 
more accurate results when we use non-interacting cluster wave functions to determine h > ( vall , 
the radius where the wave function vanishes. We also tried adding an attractive potential 
well in addition to the spherical hard wall boundary in order to be able to continuously vary 
the scattering energy. However, we found that this potential well distorts the wave function 
in the asymptotic region and results in larger errors of the phase shifts. 


We then considered the dimer-fermion system in three spatial dimensions, which also corre¬ 
sponds to neutron-deuteron scattering in the spin-quartet channel. For this calculation we 
use adiabatic projection and the spherical wall method with R' waU determined from the non- 
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interacting wave function. We find that this method gives good results for the phase shifts 
which agree well with the exact lattice phase shifts determined from Lanczos calculations of 
the spectrum of the microscopic Hamiltonian H. 

In this paper we have shown that very-high-precision energy calculations are not needed for 
determining phase shifts, and one can instead use spherical wall boundaries to measure phase 
shifts directly from scattering cluster wave functions. The methods we have presented here 
are specially designed to be immediately useful for large-scale calculations of cluster-cluster 
scattering using lattice Monte Carlo. Since the writing of the original draft of this work, the 
methods described here have been combined with lattice Monte Carlo simulations to produce 
the first ab initio calculation of alpha-alpha scattering [45]. The adiabatic projection method 
is used to reduce an eight-body system of nucleons to a system of two alpha particles. There 
has also been significant improvements on the extraction of lattice phase shifts using the 
spherical wall method with auxiliary potentials such as attractive wells and complex-valued 
potentials for particles with spin and partial-wave mixing |3K] . 

We are confident that more applications are possible which combine the adiabatic projec¬ 
tion method and Monte Carlo methods for scattering and reactions over a diverse range of 
few- and many-body systems. We note, for example, recent developments using impurity 
lattice Monte Carlo [23 S3, which opens the possibility of ab initio calculations of impurity 
quasiparticle scattering in quantum many-body systems. 
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